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ABSTRACT 

Based largely on energy budget considerations and the observed cosmic-ray (CR) ionic composition, super- 
nova remnant (SNR) blast waves are the most likely sources of CR ions with energies at least up to the "knee" 
near 10^^ eV. Shocks in young shell-type TeV-bright SNRs are surely producing TeV particles, but the emis- 
sion could be dominated by ions producing Tr'^-decay emission or electrons producing inverse-Compton (IC) 
gamma-rays. Unambiguously identifying the GeV-TeV emission process in a particular SNR will not only help 
pin down the origin of CRs, it will add significantly to our understanding of the diffusive shock acceleration 
(DSA) mechanism and improve our understanding of supernovae and the impact SNRs have on the circum- 
stellar medium. In this study, we investigate the Vela Jr SNR, an example of TeV-bright non-thermal SNRs. 
We perform hydrodynamic simulations coupled with non-linear DSA and non-equilibrium ionization near the 
forward shock (FS) to confront currently available multi- wavelength data. We find, with an analysis similar to 
that used earlier for SNR RX J 17 13.7-3946, that self-consistently modeling the thermal X-ray line emission 
with the non-thermal continuum in our one-dimensional model strongly constrains the fitting parameters, and 
this leads convincingly to a leptonic origin for the GeV-TeV emission for Vela Jr This conclusion is further 
supported by applying additional constraints from observation, including the radial brightness profiles of the 
SNR shell in TeV gamma-rays, and the spatial variation of the X-ray synchrotron spectral index. We will 
discuss implications of our models on future observations by the next-generation telescopes. 

Subject headings: acceleration of particles, shock waves, ISM: cosmic rays, ISM; supernova remnants 



1. INTRODUCTION 

Young shell-type supernova remnants (SNRs) with strong 
non-thermal emission, such as RX J1713. 7-3946, HESS 
J1731-347 and Vela Jr, have been in the spotlight in recent 
years due to their presumed close relation to the origin of 
Galactic cosmic rays (CRs). These remnants show high lu- 
minosities in TeV gamma-rays and synchrotron dominated X- 
rays with no sign of thermal emission from the ejecta or shell. 
Vela JrQ is a typical example of a TeV-bright SNR where 
emission in the X-ray band is strongly non-thermal. It lies 
along the line-of-sight of the crowded Vela complex, compris- 
ing the Vela SNR, the Vela-X pulsar wind nebula (PWN), the 
Vela pulsar, the Pencil Nebula, and the Puppis A SNR. While 
Vela Jr was first discover ed by the ROSA T satellite at X-ray 
energies above ^ 1.3 keV (I AschenbachI 1 9 98). it is coincident 
with the southeastern part of the more diffuse Vela SNR and 
is totally obscured by the latter in the soft X-ray band. The 
age of and distance to this remnant have been debated with 
the current most-agreed-on v alues being d.sNR 750 pc and 
tage - 1700 - 4300 yr (e.g.. lSlane etll]|200ll: iKatsuda et alj 
120081) . The progenitor has been suggested to be the core- 
collapse of a massive star due to the discovery of a central 
compact object AX J0 851. 9-4617.4 near the center of the 
SNR (Slane et al. 1200 Ih with a consistent age and distance 
(iKargaltsev et alJ 20021) although no pulsation has yet been 
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detected. 

Detection of the 1.157 MeV radioactive decay line of the 
short-lived "^^Ti alo ng the same line-of-sight has been claimed 
jlvudin et alj|1998i) , suggestive of an extremely short distance 
of about 200 pc, and a young age of around 700 yr How- 
ever, the line detection s uffers from low statistic al signifi- 
cance at the 2 — 4(7 level (ISchonfelder et al.ll2000h . and from 
the non-detection of the associated Ca(Sc) Ka line in X-rays 
by the Adva nced Satellite for Cosmology and Astrophysics 
(ASCA ) GIS jSlane et al.ll2001l) and Suzaku XIS jHiraga et alj 
120091) imaging spectrometers. Moreover, using the distance 
and age deduced dir ectly from proper-motion measurements 
(iKatsuda et al.l 120081) . it can be shown that the required ini- 
tial mass of ^''^Ti has to be ~ 1 Mq, orders of magnitude 
larger than that expected by supernova nucleosynthesis mod- 
els. Hence, we consider the association of the "'^Ti flux with 
Vela Jr unlikely. 

Sharp synchrotron filamentary structures in X-rays have 
been resolved by C handra in the bright northwestern shell 
(iBamba et al.l 120051) . An explanation of the small effective 
width of these filaments by fast synchrotron cooling of CR 
electrons req uires a local rnagnet ic field strength of about 
100 yuG (e.g.. lBerezhko et al.l2009h . probably indicating mag- 
netic field amplifi cation (MFA) induced by efficient CR ion 
acceleratio n (e.g.. lBellll2004t IVladimirov. Bykov. & EUisonl 
|2008; Bvko v Osipov. &Ellisonll2qill). . . 

The radio remnant (ICombi et al.l 1 19991: IStupar et al.ll2005l) 
is relatively faint (~ 30 — 50 Jy at 1 GHz) with a flat spec- 
tral index ^ 0.3. The morphology is ring-like and shows a 
good overall correlation with the X-ray remnant, suggesting a 
common synchrotron origin for the radio and X-ray emission. 
Spatially extended emission in the GeV and TeV bands has 
also been discovered by Fermi L AT (Tanaka et al. 2011) and 
H.E.S.S. jAharonian et al.ll2007l) along the line-of-sight. The 
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gamma-ray spectra can be described by a single power-law 
with index ^ 1.85 and a cutoff energy at a few TeV. As we 
will discuss in more detail below. Vela Jr shares some strik- 
ing similarities with another well-studied non-thermal SNR 
RX J1713. 7-3946 in both spectral and morphological charac- 
teristics. We believe these similarities are important for un- 
derstanding the physical nature o f Vela Jr 

Recently, iTanaka et alj (1201 lb have fitted the broadband 
spectra of Vela Jr using both hadronic (7r°-decay dominated) 
and leptonic (IC dominated) scenarios. They assumed sin- 
gle power-law spectra with exponential cut-offs for the un- 
derlying proton and electron distributions. To obtain the best 
fits, they varied the spectral indices of the particles and the 
environmental parameters including the downstream (DS) B- 
field strength and gas density. They concluded that reasonable 
spectral fits can be obtained under both scenarios by invoking 
different parameter choices. In some cases, however, the in- 
terpretation of their parameters can be non-trivial, such as an 
unrealistically large energy in the accelerated proton popu- 
lation 10^^ erg) in the hadronic model for a low-density 
(n ^ 0.01 cm^^) ambient medium, and a small downstream 
B-field (i?2 ^ 10 fiG) required in the leptonic case which 
seems to contradict with the sharpness of the X-ray filaments 
found in the brighter NW rim (but see our arguments for such 
a possibility in Section|5]). 

While DSA models predict the acceleration of both elec- 
trons and ions, determining the total energy appearing in rel- 
ativistic particles depends crucially on the correct interpre- 
tation of which particles are responsible for any observed 
gamma-ray emission. In addition, models for which the 
gamma-ray emission is produced predominantly by hadrons 
necessarily produce observable signatures in the luminosity 
and ionization state of the thermal X-ray emission. We be- 
lieve that disentangling the emission origin requires a more 
self-consistent SNR model, as we have developed and used 
to model SNR RX J1713.7-3946 (e.g.. lElUson etaD lIoTl 
iLee et al] 120121 and references therein). In our CR-hydro- 
NEI model, the thermal X-ray emission is determined self- 
consistently using the dynamical information of the forward 
shock. Furthermore, the SNR dynamics is calculated self- 
consistently with DSA including the feedback effects from 
NL-DSA on the SNR evolution. We include spatial profiles of 
the multi-wavelength emission which provide more stringent 
constraints on the parameters and lead to a better differentia- 
tion between the 7r"-decay and IC-dominated cases. 

Within the approximations of our one-dimensional model, 
we find that leptonic models clearly fit the broadband obser- 
vations better than tt*^ -decay ones in a fashion strikingly sim- 
ilar to broadband models of SNR RX J1713.7-39460 We 
emphasize that regardless of the approximations intrinsic to 
our model; our results clearly show that including a self- 
consistent calculation of the thermal X-ray emission is essen- 
tial for any broadband model of a young SNR. For Vela Jr, as 
for SNR RX J 17 13.7-3946, hadronic scenarios are excluded 
with high confidence in homogeneous models. 

This paper is organized as follows: The first section pro- 
vides a brief description of the latest version of the CR-hydro- 
NEI simulation code we use as our modeling platform; the 



^ As in our models of SNR RX J1713. 7-3946, we only calculate emission 
from the shocked region between the contact discontinuity (CD) and FS and 
neglect emission from the reverse shock. Including additional thermal X-ray 
emission from the shocked ejecta would strengthen the case for a leptonic 
dominated model for the 7-ray emission. 



second section describes our models for Vela Jr and the vari- 
ous constrains from observation data; and in the last part, we 
discuss the implication of our models on our understanding of 
the environment of Vela Jr, the acceleration of CR particles at 
the blast wave, and the possible production mechanism(s) of 
the observed broadband emission from radio to the TeV band. 
Some concluding remarks follow. 

2. THE CR-HYDRO-NEI SIMULATION CODE 

When it comes to modeling emission from young shell-type 
SNRs with strong forward shocks, it is very important to take 
into account the non-linear aspects of DSA and its coupling to 
the hydrodynamics and plasma conditions in the shocked gas. 
The high Mach number shocks in young SNRs are expected 
to be efficient particle accelerators and with efficient DSA, 
the CR production, shock structure, shocked gas temperature, 
and magnetic fields are coupled and feedback effects cannot 
be ignored. Here we employ a CR-hydro-NEI code which can 
effectively model young, shell-type SNRs like Vela Jr with 
these non-linear effects accounted for self-consistently. In this 
section, we briefly describe a generalized version of the CR- 
hydro-NEI code that includes several updates and additional 
functio ns not included in the description given in iLee et al.1 
(l20T2h . 

2.1. Coupled SNR hydrodynamics, NL-DSA, andNEI 

The development of the CR-hydro-NEI code used here has 
been described in detail in a number of recent papers (e.g., 
.Ellison et al . 2007; Patnaude et al. 2009, 2010; EllisonSaD 
I2OIOI 2012h. with th e most recent "generalized" version given 
in iLee et alJ (120121) . Briefly, the SNR hydrodynamics are 
modeled wi th a one-dimensional hyd ro simulation based on 
VH-I (e.g.. iBlondin & Ellison|[200lh . The SNR simulation 
provides the evolving shock speed, sonic and Alfven Mach 
numbers, and other quantities necessary for the shock accel- 
eration calculation. The NL-DSA calculation is done with a 
semi-analytic solution based largely on the work of R Blasi 
and co-workers (e.g ., iBlasi et al.|[2005t ICaprioU et aP 120091) 



(see lLee et al.l2012l for a full list of references). With suitable 
parameters, the semi-analytic calculation provides the shock 
compression ratio, amplified magnetic field, and full proton 
and electron spectra which are used to calculate the broad- 
band continuum radiation from synchrotron, bremsstrahlung, 
IC, and 7r°-decay emission. 

The SNR evolution is coupled to the CR production by 
modifying the standard hydrodynamic equations through a 
change in the equation of state from the influence of rel- 
ativistic CRs, energy loss from escaping CRs, and mag- 
netic pressure^ The changes produced in the hydrodynam- 
ics by NL-DSA also modify the evolution and ionization state 
of the shocked thermal plasma and these changes are self- 
consistently included in our non-equilibrium ionization (NEI) 
calculation of the thermal X-ray emission. 

We thus obtain a model of an evolving SNR where the non- 
thermal continuum and thermal line emission are calculated 
self-consistently. Treating thermal and non-thermal processes 
consistently is required when CR production is efficient be- 
cause the production of relativistic CRs can strongly influ- 
ence the density and temperatur e of the shock-heated ther- 
mal plasma (e.g., lEllisonI I2OOOI) . 



In such a self-consistent 



* We note that only a "scalar" magnetic pressure is used in our 1 D hydro 
simulation. 
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calculation, it is not possible, for instance, to adjust param- 
eters so that TT^-decay matches GeV-TeV observations with- 
out modifying the fit to thermal X-ray line emission. As we 
have demonstrated i n our fits to SNR RX J1713.7-3946 (e.g., 
Ellison et al.l 120 12t iLee et all 120121) . this coupling strongly 
constrains broadband emission models. 

2.2. Precursor CR populations and emission 

To determine the CR proton distribution function f'^'^{x,p) 
at a certain position x in the precursor u pstream of t he for ward 
shock, we follow the recipe described in ILee et ^11(120121) [see 
equation (17)] and references therein. As for the electrons , we 
implicitly assume that they do not affect the shock structure 
and dynamics due to their low energy density compared to the 
protons, which we have checked to be the case a posteriori for 
all models presented here. We also assume that radiative cool- 
ing of the electrons occurs mainly in the downstream region 
where the B-field is highest and the electrons spend most of 
their time scattering per shock-crossing. As a result, this al- 
lows us to calculate the electron distribution in the precursor 
easily using the precursor profile and the distribution function 
at the subshock in the same way as for the protons after the 
DSA solution has been obtained at every time step. 

At each position x in the precursor, we read off the lo- 
cal environment variables such as the gas density, the B-field 
strength and the temperature as provided by the DSA solution 
to calculate the local photon emissivity. The total emission 
spectrum from the CR precursor is thus obtained by an inte- 
gration over the precursor volume. 

2.3. Secondary particle production and emission in the SNR 

shell 

Accelerated protons and heavier ions will interact with the 
background gas and produce charged pions through p-p in- 
teractions. These pions further decay into secondary parti- 
cles such as e+, e~ and v's. The charged secondaries can 
accumulate in the SNR shell and radiate as the shock propa- 
gates into the upstream medium. We have implemented this 
process in the CR-hydro-NEI code by following the contin- 
uous production of the secondaries in the post-shock region, 
as well as the adiabatic and radiative losses of the e+-e^. 
We only consider secondary production and radiation in the 
post-shock region because the density and magnetic field are 
considerably higher than in the shock precursor The relevant 
inclusive secondary cross-sections are obtained using the pa- 
rameterized model by Kamae et al. (2006t), as we do for the 
gamma-rays. Contribution from heavy ions like p-He is taken 
into account via an effective enhancement factor. These sec- 
ondaries contribute to the photon emissivity via synchrotron, 
IC and bremsstrahlung and is estimated at the end of a simu- 
lation. 

We note that there is a possibility that these secondary e+- 
e^ can be injected into DSA and be 'reaccelerated' in the 
same manner as the primary electrons. This is especially true 
for those being generated in the CR precursor However, since 
most of the secondary particles are produced and accumulated 
in the DS region, and the probability of these populations trav- 
eling back to the shock from their production sites in order 
to be injected into acceleration is relatively low, we will ig- 
nore the process of reacceleration. We will show that, whilst 
the secondary particles do not contribute significantly to the 
volume-integrated broadband photon spectrum for all models 
presented in this work, the production of these particles in the 



shocked gas and their accumulation can have important ef- 
fects on the spatial variation of synchrotron emission in some 
specific X-ray bands for a hadronic model. This result will be 
discussed in detail in Section |43] 

2.4. Multi-wavelength emission profiles 

Using the spatial information from the CR-hydro-NEI cods, 
it is relatively straightforward to obtain model predictions for 
multi-wavelength emission profiles. Line-of-sight projection 
effects are calculated assuming spherical symmetry. The pro- 
jected radial (i.e., outward from the SNR center) profiles, in- 
tegrated over a set of chosen wavebands, are convolved with 
gaussian kernels to mimic the point-spread-functions of vari- 
ous instruments and facilitate direct comparison with data. 

3. MODELS 

The radio-to-TeV spectral and morphological features of 
Vela Jr largely resemble SNR RX J1713.7-3946. The ob- 
served similarities between Vela Jr and SNR RX J1713.7- 
3946 include: (i) relatively dim radio emission with a flat 
spectral index; (ii) radio and X-ray emission which are pre- 
dominately synchrotron in origin with ring-like morpholo- 
gies, and with no evidence of thermal emission from either 
the shell or the ejecta (central) region; (iii) luminous in TeV 
energies, but relatively faint in the GeV band when compared 
to other gamma-ray SNRs, especially the middle-aged rem- 
nants h ke W51C, W44, and IC 443 (e.g.. lAbdo etallllOOgl 
I2010allbl) ; and (iv) TeV emission is spatially resolved to a 
shell-like structure and correlates well with the X-ray and ra- 
dio images. As might be expected, the observational similar- 
ities between Vela Jr. and SNR RX J1713.7-3946 lead natu- 
rally to a similar conclusion from our CR-hydro-NEI model 
that IC dominates the GeV- TeV emission process. 

3.1. Input parameters 

We have constructed two representative models: a 
"hadronic" model where the GeV- TeV emission is 7r"-decay 
dominated, and a "leptonic" model where it is IC dominated. 
For both models we assume: (i) the associated supernova (SN) 
had a kinetic energy of _EgN = 10^^ erg; (ii) the age of the 
remnant is 2500 yr to match the observed size of the SNR at 
the distances assumed in Table [T] (iii) the injection parameter 
which determines the acceleration efficiency of DSA is set at 
Xinj = Pinj/pth + ^2 / c = 3.6, coiTcsponding to a fraction 
77 ^ 2 — 3 X 10^^ of thermal particles being injected in DSA 
at any time (note that our Xinj is C Blasi et al. 2005); (iv) 
we have conservatively assumed that the electron temperature 
equilibrates with the proton temperature through Coulomb 
collisions unless otherwise mentioned; (v) the spectrum of 
the background photon fields that are up-scattered by the rel- 
ativistic electrons t o gamma-ray energies (IC) is taken from 
iPorter et all (120061) at the position of the SNR, which con- 
tains the cosmic microwave background (CMB) and the local 
interstellar radiation fields in the infrared and optical bands; 
(vi) the upstream gas temperature is fixed at Tq = 10"' K re- 
gardless of the gas density gradient; (vii) a number fraction of 
helium /ho = 0.0977 is specified for the plasma, which con- 
tributes additionally to the various photon emission. However, 
other than this scaling, we ignore the acceleration of heavy 
ions in this paper; (viii) we assume that the ambient B-field 
is oriented quasi-parallel to the shock normal over the whole 
surface, that is, we assume that the magnetic field geometry 
is unimportant as would be the case for efficient DSA where 
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the self-generated magnetic turbulence mediating DSA is to- 
tally tangled in the CR precursor and downstream from the 
shock, i.e., the Bohm limit for C R diffusion is obt ained. For 
full details of the formulation, see lLee et al.l (120121) . 

An important difference between our hadronic and leptonic 
models is that the hadronic model evolves in a uniform up- 
stream medium, typical of a Type la SN, and the leptonic 
model evolves in a wind cavity, typical of a core-collapse SN. 
While the observational evidence favors a core-collapse origin 
for the Vela Jr SNR, we could not find a set of core-collapse 
parameters that would produce a reasonable fit to the broad- 
band data for a hadronic model radiating in a pre-SN wind. 
This is because the hadronic model needs a high ambient gas 
density and high B-field strength in order to simultaneously 
explain the radio/X-ray synchrotron emission and the gamma- 
ray flux. We thus used a uniform ISM model for our hadronic 
fit. On the other hand, a leptonic model requires a lower am- 
bient gas density to suppress the 7r°-decay gamma-rays, and 
this is natually realized by a wind bubble environment for 
a remnant sitting close to the Galactic plane. More details 
will be given in Section |4] where we will further discuss this 
point. Consistent with the different SN scenarios, we assume 
an ejecta mass Mcj = 1.4 M© for our hadronic model and 
Mcj = 3 M0 for our leptonic model. 

In Table [T] we summarize the main input parameters and 
output quantities measured at the end of the simulation for 
each model. We constrain the input parameters via a satis- 
factory fit to the observed broadband spectrum from radio to 
TeV energies, as well as a general agreement of the dynami- 
cal variables (e.g., radius of the forward shock) with current 
ob servations. The distan ce to the SNR is recently constrained 
bv iKatsuda eTaTI (I2OO8I) usmg proper-motion measurements 
over a 7-year time period. They provided the estimate equiv- 
alent to dsNR > (750 ± 205) x (i;sk/30GG km s^^) pc. We 
will adopt dsNR for our models which is consistent with this 
limit. 

4. RESULTS AND DISCUSSION 
4.1. Dynamics 

We plot the time evolution of selected important physical 
quantities for the hadronic and leptonic models in Figures [T] 
and|2] The very different undisturbed upstream environments 
give rise to large differences in the dynamical behavior of the 
two models. With the lower B-field and gas density in the 
wind cavity for the leptonic model, the following are expected 
and observed compared to the hadronic model: (1) the shock 
sweeps up material at a much lower rate, so that the shock 
speed and Mach numbers decay slower, and shock radius in- 
creases faster, with time; (2) the acceleration efficiency, and 
hence the fraction, /sn, of supernova explosion energy, E^n, 
converted into CR particles at the current age is much lower 
(/sN — 0.14 vs 0.48) due to the slower injection rate of ther- 
mal particles from the shocked gas into DSA. This is also re- 
flected in the much more moderate shock modification (i.e., 
smaller i?tot) than for the hadronic model on average. The 
high fs N fraction of the hadr onic model is consistent with the 
result of iTanaka etaH (1201 Ih : (3) the lower B-field in the DS 
region {B2) results in a longer synchrotron loss time-scale, in 
this case longer than the acceleration time-scale throughout 
the age of 2500 yr, resulting in a common Pmax for electrons 
and protons. On the contrary, the hadronic model has a much 
higher ambient B-field strength in the uniform ISM, such that 
Pmax of electrons becomes limited by the synchrotron loss 
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Fig . I . — Time evolution of output quantities for our leptonic model (solid 
lines) and hadronic model (dashed lines). From the top panel: (1) shock 
radius -Rps; (2) shock velocity Vs]^; (3) sonic Mach number A/g (black) and 
Alfven Mach number M\ (red) of the forward shock; and (4) B-field strength 
in the far upstream Bo (black), right in front of the subshock Bi (green), and 
immediately behind it B2 (red). 

(and IC loss to a lesser extent) starting from about 200 yr 
These evolutionary differences relate directly to the broad- 
band emission of the SNR, as we will discuss in the upcoming 
sections. 

4.2. Broadband Photon Spectrum 

As mentioned in Section [3j to obtain reasonable fits to the 
observed photon spectra, the two models require that the SNR 
evolves within very different surrounding environments. For 
the hadronic model, in order for the 7r°-decay photons to dom- 
inate over the IC component in the gamma-ray band, the gas 
density has to be high and the number ratio between electrons 
and protons (A'cp) has to be particularly low {Kcp ^ 10~^). 
The high density requirement argues against the SNR prop- 
agating into a pre-SN wind with rapidly decreasing density 
and we are forced to assume the blast wave runs into an inter- 
stellar medium (ISM) with uniform gas density and magnetic 
field strength. 

In constrast, the leptonic model requires that the forward 
shock propagates into a magnetized pre-SN wind region. This 
introduces a descending gradient for the gas density and B- 
field. The low shocked gas density in the downstream photon- 
emitting region suppresses the 7r°-decay and bremsstrahlung. 
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Fig. 2. — Continuation of Figure [T] where again, solid curves are our lep- 
tonic model and dashed curves are our hadronic model. From the top panel: 

(1) compression ratios at the subshock -Rsub (red) and total iJtot (black); 

(2) fraction of Esn converted into relativistic particles /sn for the total ac- 
celerated (black) and escaping (red) populations; (3) acceleration efficiency 
edsa; and (4) maximum momentum, Pniax. attained by the accelerated pro- 
tons (black) and electrons (red). 

as well as the X-ray line emission, relative to the IC compo- 
nent. 

Figure |3] shows the corresponding fits of the broadband 
spectra to currently available observations. We see that both 
models can explain the radio and GeV-to-TeV gamma-ray 
data reasonably well under the respective choices of parame- 
ters. In the X-ray band, however, the hadronic model has diffi- 
culty reproducing the X-ray photon index observed by ASCA, 
mainly due to the much faster synchrotron loss rate caused by 
the higher B-field in the hadronic model. The higher ambient 
gas density required for tt" -decay dominance in the gamma- 
ray band is also accompanied by a high flux from thermal X- 
ray lines, which exceeds the X-ray flux observed by ASCA in 
the soft X-ray energy range below a few keV, despite the adop- 
tion of a conservative model of electron heating via Coulomb 
collisions in the post-shock region. The self-consistently cal- 
culated thermal X-ray emission provides a strong additional 
model constraint. 

Another important difference between the two models is 
the overall normalization factor required for the broadband 
emission to match the observed flux. For the leptonic model 
this factor is close to unity, but the hadronic model has to be 
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Fig. 3. — Lower panel: Broadband SED fit with parameters tuned so the 
gamma-rays are dominated by hadronic emission. The solid lines show emis- 
sion spectra produced by the primary CRs in the post-shock region, with the 
exception of the yellow line which shows emission from the interaction of 
the escaped CRs with the upstream medium beyond the free escape boundary 
(FEB). The dashed lines show emission in the post-shock region produced 
by the secondary e+ and e^ originating from p-p interactions and the subse- 
quent decay of charged pions. The dash-dotted fines sfiow emission from the 
CR precursor region ahead of the forward shock. Upper panel: Same as the 
lower panel but with parameters tuned so that the gamma-rays are dominated 
by IC emission. All model spectra in the lower panel are boosted by a factor 
of~ 5 to match the observed Huxes. (Data: black - radio ( Combi et al. 1999]); 
red ba nd - ASCA GIS nAharonian et al.i2Q07l): red - Fermi LAT cTanaka et all 
I2OT1I) : green - //.g.5.5 jAharonian et al.l2007l) . 

boosted by a factor ~ 5 to explain the observed flux level. 
This factor ^ 5 is included in all of the model spectra shown 
in the lower panel of Figure |3] 

4.3. X-ray Spectrum and Its Spatial Properties 

A more in-depth comparison with the currently available 
X-ray data is valuable to differentiate SNR models. Firstly, to 
see the general agreement of the models with data, we plot the 
calculated non-thermal and thermal X-ray spectra against the 
ASCA result (e.g., Slane et al. 2001; Aharonian et al. 2007) in 
Figure |4] For the hadronic model, to show the possible range 
of variance of the thermal spectrum based on different equili- 
bration models for the electron temperature, we also show the 
calculated thermal spectrum assuming that the electron tem- 
perature equilibrates instantaneously with the proton temper- 
ature behind the FS. Coulomb equilibration and instant equi- 
libration are two extremes for electron heating and the actual 
heating is very likely to lie somewhere in between. For ei- 
ther extreme in heating, the hadronic model predicts a thermal 
X-ray flux substantially higher than what observation implies 
(especially below ~ 2 keV), in addition to its failure to repro- 
duce the non-thermal spectral index. The leptonic model, on 
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Fig. 4. — X-ray spectra for our leptonic model (upper panel) and hadronic 
model (lower panel). In both panels, the blue and black solid lines con'espond 
to the synchrotron (including secondary contribution if any) and thermal line 
emission respectively, while the dashed hnes are the thermal continuum spec- 
tra. The red band is from ASCA GIS observation. In the lower panel for 
the hadronic model, the light grey lines also show the thermal spectra cor- 
responding to a model with instantaneous equilibration between the electron 
and proton temperatures right behind the subshock, but is otherwise the same. 
In the top panel for the leptonic model, the grey line shows the total spectrum. 
As in Figure [3] the hadronic model spectra are multiplied by ~ 5 to match 
the observed broadband flux level. 

the other hand, shows a general good fit to the data without an 
over-prediction of the thermal e mission. 

In this regard, we notice that iBerezhko et al.l (l2009h have 
proposed a hadronic model in which the shock is running into 
a wind bubble region created by the progenitor star embedded 
inside a high density gas cloud. We do not consider such a 
situation in this work since, in such a high density environ- 
ment (jigas ^ 0.01 cm^'^), the resultant thermal X-ray lines 
are expected to violate the observation data even more than is 
the case for our hadronic model. 

Secondly, we consider the X-ray spectrum near the bright 
NW rim of the remnant. Vela Jr was observed with the Chan- 
dra X-ray Observatory on 2003 January 5 and 6 (ObsIDs 3446 
and 4414) using the Advanced CCD imaging Spectrometer 
(ACIS). Standard cleaning and data reduction were performed 
using CIAO version 4.4. The merged observations yielded a 
total exposure time of 73.9 ks. Spectra, extracted using the 
specextract software, were obtained from a narrow region en- 
compassing the forward shock of the SNR (solid box in Fig- 
ure Is) and from a background region (dashed box) used to 
account for projected emission from the Vela SNR as well as 
internal. Galactic, and extragalactic background emission. 

The predicted thermal and nonthermal emission from CR- 
hydro-NEI for the projected spectral region indicated was 
converted into a table model for use in xspec for both the 
hadronic and leptonic models. The tbabs model for interstel- 
lar absorption was applied, and the model was convolved with 
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Fig. 5. — Chandra image of the NW region of Vela Jr The regions used 
for source and background spectra are shown by the solid and dashed boxes 
respectively. A circle containing a faint point source falling within the source 
region was excluded in the spectral extraction. 

the instrument response function and fit to the data. To ac- 
count for possible excess foreground emission from the Vela 
SNR, a thermal nonequilibrium ionization model {nei) with a 
distinct tbabs absorption component was included. 

The best-fit results are shown in Figure |6] where the lep- 
tonic model fit is shown as a black histogram, and the 
hadronic model is shown in red. The dashed histogram corre- 
sponds to the additional tbabsxnei model. It is clear that the 
hadronic model provides a poor representation of the X-ray 
data, both at low energies, where the predicted line emission 
is not observed, and at high energies, where the model is much 
steeper than the observed spectrum. The leptonic model pro- 
vides an excellent fit to the data (xt = 1-05 for 327 degrees 
of freedom). The column density is 7ih = (4.3 ± 0.2) x 
10^^ cm~^, in good agreement with p revious measurements 
(ISlane et al.ll2"00U iPannuti et alJl20IOli and the accompanymg 
thermal component has a temperature of 0.26lQQg keV with 
a column density of (0 — 3.3) x 10^^ cm~^, consistent with 
observed properties for soft thermal emission from the Vela 
SNR. The best-fit normalization for the CR-hydro-NEI com- 
ponent is about 15% higher than the ratio of the spectral ex- 
traction region length to the entire SNR circumference, which 
is reasonable given that the rim is somewhat brighter at this 
position (see Figure|5]l. 

Thanks to the high angular resolution of Chandra, it is also 
possible to study spat ial variation of the spectral properties 
near the FS (see, e.g.. IPannuti et al.ll2010l) . Correspondingly, 
we calculate the synchrotron photon index for a series of line- 
of-sight projected regions from the models which have differ- 
ent radial distances from the FS. The results are shown in Fig- 
ure|2] Besides the absolute difference of the index values, the 
hadronic and leptonic models show very distinct trends of the 
index variation as a function of distance from the FS. For our 
leptonic model, the synchrotron spectrum softens as we look 
inward from the FS, whereas our hadronic model exhibits a 
much weaker dependence on position. 

The major factors which determine the spectral index of the 
non-thermal X-rays at a certain distance from the SNR cen- 
ter are the following: (1) Energy loss history of the primary 
accelerated electrons, which is critical to the high-energy cut- 
off of the primary synchrotron component; and (2) the rela- 
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Fig. 6. — Chandra spectrum of the NW rim compared with the spectra 
from our leptonic (black) and hadronic (red) models. The model spectra are 
folded with the instrument response function and the fits reveal an interstellar 
absorption with a column density of njj ~ 4 X 10^^ cm^. An additional 
thermal component (dashed histograin) with Te ~ 0.3 keV and nu ~ 2 X 
10^^ cm^ is added to the leptonic model to explain the background residual 
at low energies. 

tive importance of the secondary e+-e^ to the primary elec- 
trons at that position. These secondaries are produced by the 
trapped accelerated protons interacting with the shocked gas, 
which lose energy radiatively like the primaries and mean- 
while accumulate in number with time. Different from the 
primaries, which are accelerated by the shock and then ad- 
vected downstream without replenishment, the secondaries 
at each position in the shocked gas are continuously gener- 
ated by the advected high-energy protons whose spectra do 
not evolve with time significantly other than from adiabatic 
losses. The spectral shapes of the secondary and primary syn- 
chrotron components are hence typically different. For our 
leptonic model, the effect from secondaries is obviously small 
due to the very low ambient gas density in the wind bub- 
ble. As a result, the energy loss of the primaries dominates 
the spatial behavior of the synchrotron radiation. Since pri- 
mary electrons accelerated by the shock at an earlier phase 
and now residing further inward in the SNR shell have suf- 
fered from synchrotron/IC and adiabatic losses for a longer 
period of time than those just being freshly accelerated, the 
cutoffs of their spectra are found at lower energies. In addi- 
tion, due to the falling B-field strength with radius in our wind 
model, those primary electrons accelerated earlier have expe- 
rienced a higher synchrotron loss rate on average than those 
residing closer to the FS at the current age. These lead to the 
softening of the synchrotron spectrum as shown by our result. 

For the hadronic model, the secondaries cannot be ne- 
glected. In the volume integrated photon spectrum (Figure|3j, 
the primary component dominates over the secondary one 
until around 10 keV. However, by looking deeper towards 
the SNR center where the number of secondary particles is 
larger due to accumulation and the primaries have a lower 
cut-off energy from a longer period of radiative loss, the sec- 
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Fig. 7. — Line-of-sight dependence of the non-thermal X-ray spectral in- 
dex. The synchrotron indices are extracted from regions with various radial 
distances from the SNR center and in the energy range of 2 — 10 keV. The 
black thick lines are the result from our best-fi t models. Red points are spec- 
tral fitting results from Kishishita et al] 120131) using archival XMM Newton 
data in the same energy rang e. Th e thin grey lines correspond to the modified 
models described in Section 1431 below. Note that the contact discontinuity 
lies at R/Rpg = 0.81 — 0.82 for both hadronic and leptonic models. 

ondary spectrum can dominate over the primary spectrum 
at a few keV and harden the resultant total spectrum in the 
1 — 5 keV energy range. We find that this effect can com- 
pensate the energy loss of the primaries in the SNR shell and, 
for the relatively moderate ambient gas density we assume, 
lead to a weak variation of the synchrotron index with ra- 
dius after projecting our spherically symmetric model along 
the line-of-sight. With a higher ambient gas density, a spec- 
tral hardening can be expected as we look inward from the 
shock. Usin^ Chan dra and XMM Newton dat a respectively, 
iPannuti et alT ( l20IOl) and lKishishita et alj ( 120 1 3[) have recently 
discovered a trend of softening for the non-thermal spectrum 
as the line-of-sight is moved inward from the FS towards the 
SNR center This trend provides further support for our lep- 
tonic scenario. For a direct comparis on, we overlaid i n Fig- 
ure|2]the fitted photon indices given bv lKishishita et al.l (120 1 3h 
and found a reasonable agreement with our result. Addition- 
all y, the phenom enological spectral evolution model invoked 
by lKishishita et al . (2013) shows that the typical downstream 
B-field should be around a few /iG in order to explain the 
observed softening trend, which is highly consistent with the 
field strength obtained in our leptonic model. 

4.4. Gamma-ray Morphology 

The current generation of ground-based Imaging Atmo- 
spheric Cherenkov Telescopes (lACTs) possesses sufficient 
angular resolution to discern the spatial structures of SNRs 
with large angular sizes in sub-TeV to mu hiple TeV ener- 
gies. Recent HESS observations of Vela Jr (lAharonian et al.l 
l2007h have measured the radial surface brightness profiles of 
its shell-like emission from 300 GeV to 20 TeV with a res- 
olution of about 0.06°. We compare our model profiles with 
the HESS measurements for the brightest northern part in Fig- 
ure |8] The profiles are extracted from the model gamma-ray 
image after convolution with the PSF of the observation. The 
results show that both models predict gamma-ray emission 
profiles reasonably compatible with the data, and unlike the 
X-ray result, closely resemble each other after PSF smooth- 
ing. 
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Fig. 8. — Radial profile of gamma-ray emission (300 GeV - 20 TeV) com- 
pared with HESS measurements of the northern shell. The black thick curves 
are model profiles smoothed by a Gaussian kernel with cr = 0.06° corre- 
sponding roughly to the PSF of the HESS observation. The red thin curves 
are the model profiles expected under the advertised best angular resolution 
of CTA (0.02°). 



The reason behind this relative similarity in the gamma-ray 
band can be explained by the following: X-ray synchrotron 
emission produced by relativistic electrons is highly sensitive 
to the radiative loss rate. Electrons experience much faster 
energy loss in our hadronic model than the leptonic model, 
as mentioned above, which results into a sharp X-ray profile 
close to the FS for the hadronic case, and a more spread out 
profile behind the FS for the leptonic case. In the gamma-ray 
band, however, radiative loss plays a much less important role. 
In the hadronic model, the gamma-ray photons mainly origi- 
nate from CR protons where losses can be ignored, while in 
the leptonic model the electrons responsible for the IC emis- 
sion are only subject to minor losses due to the low magnetic 
fields. 

Nevertheless, the gamma-rays originate from totally differ- 
ent mechanisms for the hadronic and leptonic cases so differ- 
ences are predicted for more precise measurements. Future 
ground-based gamma-ray observatories, such as the upcom- 
ing Cherenkov Telescope Array (CTA), will be able to differ- 
entiate the models with unprecedented spatial resolving power 
and sensitivity. For illustration, we include in Figure |8] our 
model profiles using the advertised best angular resolution of 
CTA (about 0.02°). The red curves imply that future gamma- 
ray observations may be able to differentiate between the pre- 
dicted radial profiles of Vela Jr 

4.5. Possible Thermal X-ray Line Detection by Astro-H 

TeV-bright SNRs Hke Vela Jr and RX J1713.7-3946 have 
shown no sign of strong thermal emission with current X-ray 
observations. However, with the advent of future instruments 
of superior spectral resolution, it is very possible that thermal 
X-ray lines will be discerned in the non-thermal dominated 
spectrum. As a first look, we plot the calculated X-ray spec- 
trum in Figure |9] for our leptonic model in the energy range 
of 0.3 to 12 keV, the designed energy coverage of the Soft 
X-ray Spectrometer (SXS) onboard the next-generation X-ray 
space observatory Astro-H. The emission lines are filtered by 
a Gaussian kernel with a FWHM of 7 eV to mimic the target 
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Fig. 9. — Thermal emission fines and non-thermal spectrum from 0.3 to 
12 keV from the leptonic model. The model spectrum is filtered with a Gaus- 
sian kernel with AE (FWHM) = 7 eV consistent over the whole energy range 
to imitate the (minimum) target energy resolution of the SXS instrument on- 
board the next-generation X-ray observatory Astro-H. A thick cyan line is 
overlaid to show the total spectrum under a fixed resolution of 60 eV for 
comparison. 

baseline spectral resolution of SXSQ 

For comparison, we show the leptonic model spectrum with 
a resolution of 60 eV (a typical ballpark value for currently 
operating X-ray observatories). We see that using the resolv- 
ing power of SXS, the stronger sub-keV lines can be recog- 
nized very visibly on top of the synchrotron spectrum, and it 
is very possible that they will be detected by Astro-H. Detec- 
tion of such lines will be important to further constrain our 
model parameters such as the local gas compositions and the 
ion temperatures and densities. In a follow-up paper in prepa- 
ration, we will perform a more detailed study of the X-ray 
spectrum by performing formal Astro-H SXS spectral simu- 
lations based on our models, and take into account thermal 
Doppler broadening of the emission linesj^ 

4.6. Sensitivity of Results on Model Parameters 

We emphasize that while small variations in parameters in 
either of our "best-fit" models do not strongly modify the re- 
sults for that model, the parameter sets for the two scenarios 
are well separated and it is not possible to go smoothly from 
one set to the other while maintaining a satisfactory fit to the 
observations. 

To illustration this point and see how our models respond 
to parameter changes, we create modified versions for each 
of our models in which we purposely tune a set of impor- 
tant parameters so that it is positioned closer to the parameter 
space of the competing best-fit model. Two of the most im- 
portant parameters that distinguish our leptonic and hadronic 
scenarios include the electron-to-proton number ratio at rela- 
tivistic energies A'op, and the upstream gas density uq (or the 
mass loss rate dM /At and wind speed V"„ind when a pre-SN 
wind is present). For the leptonic case with a massive star 
as progenitor and a wind cavity, by recalling that the IC flux 

' We assume a pure spherical expansion of the SNR shell, such that near 
the FS (i.e., the SNR rim) the line broadening/shifting effect from expansion 
is negligible. 

For the typical post-shock temperature in our leptonic model, the broad- 
ening by thermal motion of the ions is around 1%, assuming that the temper- 
ature ratios among the ion species are mass-proportional. Therefore, roughly 
speaking, SXS with a FWHM of at most 7 eV will be able to identify thermal 
broadening of the stronger lines above approximately 0.7 keV, where there 
are a few. We will perform a more detailed calculation of the ion tempera- 
tures in our follow-up paper to better assess this possibility. 
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fic oc 7ioKcp (e.g. lEUison et al.|[2010l) . we explore an altered 
model where Kgp is lowered by a factor of 5 (i.e. 3 x 10~^) 
and the wind matter density increased by the same factor to 
roughly maintain the level of /ic0 The wind magnetization 
is fixed so as to keep the synchrotron-to-IC flux ratio /syn/ /ic 
roughly unchanged. On the other hand, we explore an altered 
model for the hadronic case where i^op is boosted up by a 
factor of 5 while uq is decreased by a factor of 1 .7, so that the 
model is now leaning toward the leptonic parameter space. 
Note that since the 7r"-decay flux Jt^o oc n^, this will decrease 
f„a and increase /ic by roughly the same factor For both 
cases, we also adjust the assumed age of the remnant, within 
the allowable range of 1700 — 4000 yr implied by observa- 
tions, from the default age of 2500 yr in order to achieve the 
same SNR angular size as the original model. Other than what 
described above, all other parameters are kept fixed. 

The effect on the total broadband spectrum is shown in Fig- 
ure[TO] We can see that the modified leptonic model produces 
a non-thermal spectrum with an acceptable fit to the data ex- 
cept the highest energy TeV points and a X-ray spectral index 
a little too soft. But most importantly, it already starts to face 
the difficulty of over-predicting the thermal emission in the X- 
ray band due to the enhanced ambient gas density, which is the 
same problem encountered by our hadronic model. Another 
surfacing problem not shown in this plot is that the spatial 
dependence of the synchrotron index now fails to explain the 
XMM Newton data by showing a softening trend too strong 
towards the SNR center (see Figure |7). For the modified 
hadronic model, although the thermal X-ray reduces back to a 
flux level more compatible with data than the original model, 
the broadband spectral fit becomes unacceptable by failing to 
explain the non-thermal X-ray flux level simultaneously with 
the gamma-rays, due to the enhanced leptonic contributions 
which worsen the fit. To achieve a good fit again, we find 
that a much lower B-field must accompany the increased A'cp 
and decreased no, which is simply what being realized in the 
wind model we are employing for the leptonic model. The 
conclusion from this analysis is hence that it is impossible to 
maintain a good fit to data by continuously tuning some influ- 
ential parameters from the leptonic to the hadronic case (and 
vice versa) unless a drastic change of the underlying SNR en- 
vironment is involved, i.e. a rarified wind cavity against a uni- 
form ISM model for the ambient medium. 

5. CONCLUSIONS 

We have investigated the TeV bright SNR Vela Jr with 
a comprehensive modeling of its multi-wavelength emission 
from the radio to TeV band. Our one-dimensional, spheri- 
cally symmetric NL-DSA model indicates that the SNR orig- 
inated from the core-collapse explosion of a massive star in- 
side a wind cavity with rarified gas density and magnetic 
fields. Based on the broadband continuum emission, with the 
added constraint from the self-consistently determined ther- 
mal X-ray emission, we show that the GeV-TeV gamma-rays 
are predominately produced by IC emission from shock ac- 
celerated electrons with a much smaller hadronic contribution 
from shock accelerated ions. Despite the fact that the leptons 
radiate more efficiently, here, as in all self-consistent DSA 
models of SNRs we are aware of, far more energy is put into 
relativistic ions than electrons. In summary, these conclusions 

' Only decreasing Kcp by such a factor will not visibly change the broad- 
band spectral shape, but will require an overall normalization factor substan- 
tially larger than 1, thus excluding such an option. 
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Fig. 10. — Comparison of the broadband spectra of our original best-fit 
models (dashed lines) with their respective modified versions (solid lines) to 
show the sensitivity of our models on parameter choices. In each panel, the 
black lines show the total non-thermal spectrum, and the grey lines show the 
accompanying thermal X-ray line emission. The data points are the same as 
Figure[3] See text for details about the modified models. 

are supported by the following results: 

• In order for 7r°-decay gamma-rays from shock acceler- 
ated ions to dominate IC emission from electrons, the 
ambient pre-shock density must be above some limit. 
Despite a wide range in possible thermal equilibration 
models for electron heating and the non-equilibrium 
ionization of the shocked plasma, the density required 
for 7r"-decay dominance is accompanied by strong ther- 
mal line emission in the X-ray band which is well above 
ASCA GIS observations (see the lower panel of Fig- 
ure [3]l. This conclusion is similar to what we found 
previously for another non-thermal SNR RX J 17 13.7- 
3946. 

• Even if the inconsistency with the thermal X-ray line 
emission is ignored, the hadronic model, with the SNR 
forward shock moving into a uniform ISM, fails to 
reproduce the spectral index of the observed X-ray 
synchrotron emission, as measured by ASCA GIS and 
Chandra ACIS. While the hadronic model can pro- 
duce a reasonable match to the spectral shapes of data 
in other wavelengths, the overall normalization of the 
broadband emission requires a boosting factor of 5 
to bring the model spectra to the observed flux levels. 

• Another possible problem with the hadronic model is 
that our best fit hadronic result predicts a spatial (i.e., 
radial) variation of the X-ray synchrotron photon in- 
dex around the FS that is incompatible to the trend sug- 
gested by data from X-ray satelUtes (see Figure |7). 
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• The Vela Jr. SNR most likely originated from a core- 
collapse SN but the "best-fit" hadronic model requires 
a relatively high-density, high-magnetic field environ- 
ment more typical of a thermonuclear Type la SN. 

In contrast to the above problems with the hadronic model, 
our pre-SN wind leptonic model can self-consistently ex- 
plain the broadband continuum emission, the lack of ther- 
mal X-ray line emission, the overall normalization assuming 
rfsNR ^ 1 kpc, and the spatial variation of the non-thermal 
X-ray spectral index. 

It is important to note that while our calculations as- 
sume a uniform ambient ISM, any real SNR may have a 
highly inhomogeneous upstream medium containing clumped 
gas. Recently, radio observations of SNR RX J1713.7-3946 
(ISano et al.l 120101: iFukui et al.l 120 12h have revealed clumpy 
atomic/molecular cloud distributions surrounding the SNR 
shell, indicating a possible interaction of the shock and the ac- 
celerated protons with the gas clouds. Concentrating the mass 
in small, dense clumps can considerably enhance the gamma- 
ray flux. Ilnoue et aiT(l2012l) have investigated this possibility 
using 3-D MHD simulations of a shock running into an inho- 
mogeneous medium. They pointed out the possibility that a 
hadronic model can avoid the problem of an over-prediction 
of the thermal X-ray flux, since the denser cores of the gas that 
carry most of the mass fraction can survive the shock pass- 
ing without being ionized and hence will not contribute to the 
thermal emission, but still can interact with the accelerated 
protons and produce gamma-rays. While definitively ruling 
out hadronic scenarios in homogeneous broadband models of 
Vela Jr, our 1-D model cannot rule out hadronic dominated 
GeV-TeV emission in more complex clumpy environments. 
In fact, although not pursued in this work, it is possible that 
the gamma-rays are contributed by a two-component mixture 
of IC and Tr^'-decay photons from the interaction of dense 
gas clumps and the accelerated CR protons. Future high- 
resolution observations of the surrounding interstellar gas dis- 
tribution and 3-D hydro simulations self-consistently coupled 
to DSA will provide a quantitative estimate on this possible 
7r"-decay component from clumps. 

To explain the broadband spectrum of the SNR as a whole, 
our leptonic model requires a low B-field a few fiG) on 
average in the SNR shell, which is one of the reasons why 
the presence of a wind cavity is necessary to provide a weak 
B-field in the upstream medium. The fact that magnetic field 
amplification in the shock precursor is taken into account in 
our models only strengthens this need. This low downstream 
B-field is highly consistent with those inf erred from mea- 
surem ents of the X-ray/TeV brightness ratio (lAharonian et alj 
l2007h and more recentl y the radial variation of X-ray syn- 
chrotr on spectral indices (iPannuti et al.l2010tlKishishita et alj 
120131) . However, as briefly mentioned earlier, it disagrees with 
some previous claims of large post-shock B-fields ('^ 0.1 mG) 
based on the discovery of sharp X-ray filaments in the NW 
shell. The argument for such large field strength is built on the 
assumption that the narrow widths of such filaments are real- 
ized by a short synchrotron loss time-scale for the electrons. 
Our broadband analysis and supporting observational data es- 
sentially show that an alternative explanation for the sharp 
filamentary structures without the need to invoke large B- 
field strengths may be nece ssary (see e.g.. lBvkov et al.luOOSi : 
lEUison & Vladimirovl2008l and references therein). We leave 
this important point as an open question to be answered by fu- 
ture works. 



A fundamental question for DSA in SNRs is the fraction 
of SN explosion energy, /sn, channeled into CRs. Our two 
models show that this fraction can be very different, at a 
given age, depending on whether the gamma-ray production 
mechanism is assumed to be 7r° -decay or IC. For Vela jr, 
at tsNR = 2500 yr, /sn — 0.14 for our leptonic model 
and /sN — 0.48 for our hadronic model. In both cases, of 
course, CR protons contain the large majority of the CR en- 
ergy. The ratio of CR electron to CR proton energy density is 
given by /\cp and Kcp — 0.015 for our leptonic model and 
A'op ~ 1.5x10^4 for our hadronic model. 

The determination of the gamma-ray origin, therefore, re- 
lates directly to the long postulated but unproved link between 
SNRs and Galactic CRs, that is, whether or not the Galac- 
tic ensemble of SNRs can meet the energy budget of the CR 
spectrum. Because the 7r°-decay and IC mechanisms are so 
different, a mis-interpretation of the multi-wavelength data in- 
cluding the gamma-rays can yield a very wrong prediction for 
the contribution of SNRs to the Galactic CRs. Of course, the 
CR energy budget depends ultimately on /sn over the full 
age of the remnant so it will be interesting to see how /sn 
evolves from very young ages up to the radiative phase of a 
SNR. In particular, the recent discovery of gamma-ray bright 
middle-aged SNRs by Fermi LAT may imply active particle 
acceleration at radiative shocks running into molecular clouds 
despite their slow speeds. 

Another important point is that the evolution of a SNR in 
a changing ambient environment may shift the dominating 
gamma-ray production mechanism with age. It is hence crit- 
ical to understand observations from SNRs of different ages 
and try to link them together into a self-consistent evolution- 
ary picture. We will pursue this line of work in a series of 
follow-up papers. 

Finally, we discussed the implications of our models for 
Vela Jr on future observations by the next-generation of ob- 
servatories, including CTA and Astro-H. We showed that CTA 
may be able to distinguish broadband emission models by 
measuring TeV brightness profiles with its unprecedented an- 
gular resolving power This is especially valuable to distin- 
guish leptonic and hadronic models for Vela Jr, as we have 
illustrated in Figure|8] In the X-ray band, the soon to be com- 
missioned Astro-H space telescope possesses an extremely 
powerful calorimeter for high-resolution spectral studies, and 
we have postulated that there is a high possibility that SXS 
will be able to detect the more prominent thermal lines (see 
Figure |9]l from non-thermal dominated SNRs such as Vela 
Jr and RX J1713. 7-3946. We expect future multi-wavelength 
observations will put extremely precise constraints on our 
model parameters, allowing us to firmly pin down the origin 
of the high-energy emission and to gain further insight on the 
particle acceleration mechanism in young SNR shocks. 
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TABLE 1 
Model Summary * 





1— TciHvonir* 


T f*r\l"onir* ^ 


IVCllltlllvft 


Input paiameters 








dsNR [kpc] 


0.74 


0.88 


Distance to SNR 


no [cm~'^] 


0.033 


(0.002) 


Far upstream gas density 


Bo [fiG] 


0.5 


(0.14) 


Far upstream B-field 


dM/dt [10-^ Mgyr-i] 





7.5 


Pre-SN mass loss rate 


V„in(i [km s~^] 





50 


Pre-SN wind speed 







0.02 


Wind magnetization 


K^p 


1.5 X 10-4 


0.015 


(e-/p) ratio at pmax of e" 


Ocut 


0.75 


0.50 


Exp. cut-off index 


^FEB 


0.15 


0.12 


Width of FEB over ijps 


J all 


0.10 


1.00 


x-dependance of v ^ 


Output quantities 








-Rfs [pc] 


12.7 


15.2 


Forward shock radius 


RcD [pc] 


10.3 


12.5 


Contact discontinuity radius 


Vps [1cm s-i] 


2130 


4700 


Forward shock velocity 


Pmax (p) [TeV/c] 


26.7 


5.2 


Proton max. momentum 


Pmax (e~) [TeV/c] 


13.3 


5.2 


Electron max. momentum 


Rtot 


9.30 


4.69 


Total compression ratio 


Rsuh 


3.69 


3.99 


Subshock compression ratio 


B2 [mG] 


34.1 


4.8 


Downstream B-field 


Ta [10« K] 


0.16 


3.62 


Downstream temperature 




0.84 


0.36 


Total accel efficiency 




0.34 


0.12 


Escaped accel efficiency 


-Ecr/-Esn (/sn) 


0.48 


0.14 


SN energy converted to CR 



* For all models listed in this Table, Bohm diffusion is assumed in the DS and precursor region, and MFA is included. These 
models all have Bsn = 10^^ erg. To = 10'* K, and /hc = 0.0977. 
t Values in parentheses are taken at t = tago (= 2500 yr). 

" Best-fit model under the 'hadronic-dominated' scenario with a uniform upstream medium. 
Best-fit model under the 'leptonic-dominated' scenario with a pre-SN magnetized wind. 
See e.g. equation (24) in Lee et al. 1 2012) for definition. 
See e.g. equation (1) an d description in Ellison & Vladimirovl 1200 8ft . 
See e.g. equation (14) in lLee et al] 120121) for definition. 



